SARS-CoV-2 mutations: the biological trackway towards viral fitness.

The outbreak of pneumonia-like respiratory disorder at China and its rapid transmission world-wide resulted in public health emergency, which brought lineage B betacoronaviridae SARS-CoV-2 (severe acute respiratory syndrome coronavirus 2) into spotlight. The fairly high mutation rate, frequent recombination and interspecies transmission in betacoronaviridae are largely responsible for their temporal changes in infectivity and virulence. Investigation of global SARS-CoV-2 genotypes revealed considerable mutations in structural, non-structural, accessory proteins as well as untranslated regions. Among the various types of mutations, single-nucleotide substitutions are the predominant ones. In addition, insertion, deletion and frame-shift mutations are also reported, albeit at a lower frequency. Among the structural proteins, spike glycoprotein and nucleocapsid phosphoprotein accumulated a larger number of mutations whereas envelope and membrane proteins are mostly conserved. Spike protein and RNA-dependent RNA polymerase variants, D614G and P323L in combination became dominant world-wide. Divergent genetic variants created serious challenge towards the development of therapeutics and vaccines. This review will consolidate mutations in different SARS-CoV-2 proteins and their implications on viral fitness.

Comparison of SARS-CoV-2 whole genome sequencing using tiled amplicon enrichment and bait hybridization.

The severe acute respiratory syndrome coronavirus 2 (SARS­CoV­2) pandemic has led to extensive virological monitoring by whole genome sequencing (WGS). Investigating the advantages and limitations of different protocols is key when conducting population-level WGS. SARS-CoV-2 positive samples with Ct values of 14-30 were run using three different protocols: the Twist Bioscience SARS­CoV­2 protocol with bait hybridization enrichment sequenced with Illumina, and two tiled amplicon enrichment protocols, ARTIC V3 and Midnight, sequenced with Illumina and Oxford Nanopore Technologies, respectively. Twist resulted in better coverage uniformity and coverage of the entire genome, but has several drawbacks: high human contamination, laborious workflow, high cost, and variation between batches. The ARTIC and Midnight protocol produced an even coverage across samples, and almost all reads were mapped to the SARS-CoV-2 reference. ARTIC and Midnight represent robust, cost-effective, and highly scalable methods that are appropriate in a clinical environment. Lineage designations were uniform across methods, representing the dominant lineages in Sweden during the period of collection. This study provides insights into methodological differences in SARS­CoV­2 sequencing and guidance in selecting suitable methods for various purposes.

Rapid assembly of SARS-CoV-2 genomes reveals attenuation of the Omicron BA.1 variant through NSP6.

Although the SARS-CoV-2 Omicron variant (BA.1) spread rapidly across the world and effectively evaded immune responses, its viral fitness in cell and animal models was reduced. The precise nature of this attenuation remains unknown as generating replication-competent viral genomes is challenging because of the length of the viral genome (~30 kb). Here, we present a plasmid-based viral genome assembly and rescue strategy (pGLUE) that constructs complete infectious viruses or noninfectious subgenomic replicons in a single ligation reaction with >80% efficiency. Fully sequenced replicons and infectious viral stocks can be generated in 1 and 3 weeks, respectively. By testing a series of naturally occurring viruses as well as Delta-Omicron chimeric replicons, we show that Omicron nonstructural protein 6 harbors critical attenuating mutations, which dampen viral RNA replication and reduce lipid droplet consumption. Thus, pGLUE overcomes remaining barriers to broadly study SARS-CoV-2 replication and reveals deficits in nonstructural protein function underlying Omicron attenuation.

The SARS-CoV-2 RNA polymerase is a viral RNA capping enzyme.

SARS-CoV-2 is a positive-sense RNA virus responsible for the Coronavirus Disease 2019 (COVID-19) pandemic, which continues to cause significant morbidity, mortality and economic strain. SARS-CoV-2 can cause severe respiratory disease and death in humans, highlighting the need for effective antiviral therapies. The RNA synthesis machinery of SARS-CoV-2 is an ideal drug target and consists of non-structural protein 12 (nsp12), which is directly responsible for RNA synthesis, and numerous co-factors involved in RNA proofreading and 5' capping of viral RNAs. The formation of the 5' 7-methylguanosine (m7G) cap structure is known to require a guanylyltransferase (GTase) as well as a 5' triphosphatase and methyltransferases; however, the mechanism of SARS-CoV-2 RNA capping remains poorly understood. Here we find that SARS-CoV-2 nsp12 is involved in viral RNA capping as a GTase, carrying out the addition of a GTP nucleotide to the 5' end of viral RNA via a 5' to 5' triphosphate linkage. We further show that the nsp12 NiRAN (nidovirus RdRp-associated nucleotidyltransferase) domain performs this reaction, and can be inhibited by remdesivir triphosphate, the active form of the antiviral drug remdesivir. These findings improve understanding of coronavirus RNA synthesis and highlight a new target for novel or repurposed antiviral drugs against SARS-CoV-2.

Current understanding of an Emerging Coronavirus using in silico approach: Severe Acute Respiratory Syndrome-Coronavirus-2 (SARS-CoV-2).

Novel coronavirus (nCoV) namely "SARS-CoV-2" is being found responsible for current PANDEMIC commenced from Wuhan (China) since December 2019 and has been described with epidemiological linkage to China in about 221 countries and territories until now. In this study we have characterized the genetic lineage of SARS-CoV-2 and report the recombination within the genus and subgenus of coronaviruses. Phylogenetic relationship of thirty nine coronaviruses belonging to its four genera and five subgenera was analyzed by using the Neighbor-joining method using MEGA 6.0. Phylogenetic trees of full length genome, various proteins (spike, envelope, membrane and nucleocapsid) nucleotide sequences were constructed separately. Putative recombination was probed via RDP4. Our analysis describes that the "SARS-CoV-2" although shows great similarity to Bat-SARS-CoVs sequences through whole genome (giving sequence similarity 89%), exhibits conflicting grouping with the Bat-SARS-like coronavirus sequences (MG772933 and MG772934). Furthermore, seven recombination events were observed in SARS-CoV-2 (NC_045512) by RDP4. But not a single recombination event fulfills the high level of certainty. Recombination mostly housed in spike protein genes than rest of the genome indicating breakpoint cluster arises beyond the 95% and 99% breakpoint density intervals. Genetic similarity levels observed among "SARS-CoV-2" and Bat-SARS-CoVs advocated that the latter did not exhibit the specific variant that cause outbreak in humans, proposing a suggestion that "SARS-CoV-2" has originated possibly from bats. These genomic features and their probable association with virus characteristics along with virulence in humans require further consideration.

SARS-CoV-2 variant identification using a genome tiling array and genotyping probes.

With over 5.5 million deaths worldwide attributed to the respiratory disease COVID-19 caused by the novel coronavirus SARS-CoV-2, it is essential that continued efforts be made to track the evolution and spread of the virus globally. The authors previously presented a rapid and cost-effective method to sequence the entire SARS-CoV-2 genome with 95% coverage and 99.9% accuracy. This method is advantageous for identifying and tracking variants in the SARS-CoV-2 genome compared with traditional short-read sequencing methods which can be time-consuming and costly. Herein, the addition of genotyping probes to a DNA chip that targets known SARS-CoV-2 variants is presented. The incorporation of genotyping probe sets along with the advent of a moving average filter improved the sequencing coverage and accuracy of the SARS-CoV-2 genome.

Throughout the COVID-19 pandemic the virus known as SARS-CoV-2 has continued to mutate and evolve. It is imperative to continue to track these mutations and where the virus has traveled to best inform healthcare practices and global strategies to combat the virus. The authors previously developed a method to investigate 95% of this viral genome with 99.9% accuracy that was more cost-effective and less time-consuming than previous methods. In this work, specific markers were added to the technology to allow tracking of mutations in the virus that have already been documented. In doing so, the accuracy and how much of the viral genome can be sequenced was improved.

Safety and Biodistribution of Nanoligomers Targeting the SARS-CoV-2 Genome for the Treatment of COVID-19.

As the world braces to enter its fourth year of the coronavirus disease 2019 (COVID-19) pandemic, the need for accessible and effective antiviral therapeutics continues to be felt globally. The recent surge of Omicron variant cases has demonstrated that vaccination and prevention alone cannot quell the spread of highly transmissible variants. A safe and nontoxic therapeutic with an adaptable design to respond to the emergence of new variants is critical for transitioning to the treatment of COVID-19 as an endemic disease. Here, we present a novel compound, called SBCoV202, that specifically and tightly binds the translation initiation site of RNA-dependent RNA polymerase within the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genome, inhibiting viral replication. SBCoV202 is a Nanoligomer, a molecule that includes peptide nucleic acid sequences capable of binding viral RNA with single-base-pair specificity to accurately target the viral genome. The compound has been shown to be safe and nontoxic in mice, with favorable biodistribution, and has shown efficacy against SARS-CoV-2 in vitro. Safety and biodistribution were assessed using three separate administration methods, namely, intranasal, intravenous, and intraperitoneal. Safety studies showed the Nanoligomer caused no outward distress, immunogenicity, or organ tissue damage, measured through observation of behavior and body weight, serum levels of cytokines, and histopathology of fixed tissue, respectively. SBCoV202 was evenly biodistributed throughout the body, with most tissues measuring Nanoligomer concentrations well above the compound KD of 3.37 nM. In addition to favorable availability to organs such as the lungs, lymph nodes, liver, and spleen, the compound circulated through the blood and was rapidly cleared through the renal and urinary systems. The favorable biodistribution and lack of immunogenicity and toxicity set Nanoligomers apart from other antisense therapies, while the adaptability of the nucleic acid sequence of Nanoligomers provides a defense against future emergence of drug resistance, making these molecules an attractive potential treatment for COVID-19.

Benchmarking machine learning robustness in Covid-19 genome sequence classification.

The rapid spread of the COVID-19 pandemic has resulted in an unprecedented amount of sequence data of the SARS-CoV-2 genome-millions of sequences and counting. This amount of data, while being orders of magnitude beyond the capacity of traditional approaches to understanding the diversity, dynamics, and evolution of viruses, is nonetheless a rich resource for machine learning (ML) approaches as alternatives for extracting such important information from these data. It is of hence utmost importance to design a framework for testing and benchmarking the robustness of these ML models. This paper makes the first effort (to our knowledge) to benchmark the robustness of ML models by simulating biological sequences with errors. In this paper, we introduce several ways to perturb SARS-CoV-2 genome sequences to mimic the error profiles of common sequencing platforms such as Illumina and PacBio. We show from experiments on a wide array of ML models that some simulation-based approaches with different perturbation budgets are more robust (and accurate) than others for specific embedding methods to certain noise simulations on the input sequences. Our benchmarking framework may assist researchers in properly assessing different ML models and help them understand the behavior of the SARS-CoV-2 virus or avoid possible future pandemics.

Genomic assessment of quarantine measures to prevent SARS-CoV-2 importation and transmission.

Mitigation of SARS-CoV-2 transmission from international travel is a priority. We evaluated the effectiveness of travellers being required to quarantine for 14-days on return to England in Summer 2020. We identified 4,207 travel-related SARS-CoV-2 cases and their contacts, and identified 827 associated SARS-CoV-2 genomes. Overall, quarantine was associated with a lower rate of contacts, and the impact of quarantine was greatest in the 16-20 age-group. 186 SARS-CoV-2 genomes were sufficiently unique to identify travel-related clusters. Fewer genomically-linked cases were observed for index cases who returned from countries with quarantine requirement compared to countries with no quarantine requirement. This difference was explained by fewer importation events per identified genome for these cases, as opposed to fewer onward contacts per case. Overall, our study demonstrates that a 14-day quarantine period reduces, but does not completely eliminate, the onward transmission of imported cases, mainly by dissuading travel to countries with a quarantine requirement.

Evidence Supporting That C-to-U RNA Editing Is the Major Force That Drives SARS-CoV-2 Evolution.

Mutations of DNA organisms are introduced by replication errors. However, SARS-CoV-2, as an RNA virus, is additionally subjected to rampant RNA editing by hosts. Both resources contributed to SARS-CoV-2 mutation and evolution, but the relative prevalence of the two origins is unknown. We performed comparative genomic analyses at intra-species (world-wide SARS-CoV-2 strains) and inter-species (SARS-CoV-2 and RaTG13 divergence) levels. We made prior predictions of the proportion of each mutation type (nucleotide substitution) under different scenarios and compared the observed versus the expected. C-to-T alteration, representing C-to-U editing, is far more abundant that all other mutation types. Derived allele frequency (DAF) as well as novel mutation rate of C-to-T are the highest in SARS-CoV-2 population, and C-T substitution dominates the divergence sites between SARS-CoV-2 and RaTG13. This is compelling evidence suggesting that C-to-U RNA editing is the major source of SARS-CoV-2 mutation. While replication errors serve as a baseline of novel mutation rate, the C-to-U editing has elevated the mutation rate for orders of magnitudes and accelerates the evolution of the virus.

Correlated substitutions reveal SARS-like coronaviruses recombine frequently with a diverse set of structured gene pools.

Quantifying SARS-like coronavirus (SL-CoV) evolution is critical to understanding the origins of SARS-CoV-2 and the molecular processes that could underlie future epidemic viruses. While genomic analyses suggest recombination was a factor in the emergence of SARS-CoV-2, few studies have quantified recombination rates among SL-CoVs. Here, we infer recombination rates of SL-CoVs from correlated substitutions in sequencing data using a coalescent model with recombination. Our computationally-efficient, non-phylogenetic method infers recombination parameters of both sampled sequences and the unsampled gene pools with which they recombine. We apply this approach to infer recombination parameters for a range of positive-sense RNA viruses. We then analyze a set of 191 SL-CoV sequences (including SARS-CoV-2) and find that ORF1ab and S genes frequently undergo recombination. We identify which SL-CoV sequence clusters have recombined with shared gene pools, and show that these pools have distinct structures and high recombination rates, with multiple recombination events occurring per synonymous substitution. We find that individual genes have recombined with different viral reservoirs. By decoupling contributions from mutation and recombination, we recover the phylogeny of non-recombined portions for many of these SL-CoVs, including the position of SARS-CoV-2 in this clonal phylogeny. Lastly, by analyzing >400,000 SARS-CoV-2 whole genome sequences, we show current diversity levels are insufficient to infer the within-population recombination rate of the virus since the pandemic began. Our work offers new methods for inferring recombination rates in RNA viruses with implications for understanding recombination in SARS-CoV-2 evolution and the structure of clonal relationships and gene pools shaping its origins.

SARS-CoV-2 detection by targeting four loci of viral genome using graphene oxide and gold nanoparticle DNA biosensor.

The current COVID-19 pandemic outbreak poses a serious threat to public health, demonstrating the critical need for the development of effective and reproducible detection tests. Since the RT-qPCR primers are highly specific and can only be designed based on the known sequence, mutation sensitivity is its limitation. Moreover, the mutations in the severe acute respiratory syndrome ß-coronavirus (SARS-CoV-2) genome led to new highly transmissible variants such as Delta and Omicron variants. In the case of mutation, RT-qPCR primers cannot recognize and attach to the target sequence. This research presents an accurate dual-platform DNA biosensor based on the colorimetric assay of gold nanoparticles and the surface-enhanced Raman scattering (SERS) technique. It simultaneously targets four different regions of the viral genome for detection of SARS-CoV-2 and its new variants prior to any sequencing. Hence, in the case of mutation in one of the target sequences, the other three probes could detect the SARS-CoV-2 genome. The method is based on visible biosensor color shift and a locally enhanced electromagnetic field and significantly amplified SERS signal due to the proximity of Sulfo-Cyanine 3 (Cy3) and AuNPs intensity peak at 1468 cm-1. The dual-platform DNA/GO/AuNP biosensor exhibits high sensitivity toward the viral genome with a LOD of 0.16 ng/µL. This is a safe point-of-care, naked-eye, equipment-free, and rapid (10 min) detection biosensor for diagnosing COVID-19 cases at home using a nasopharyngeal sample.

Global disparities in SARS-CoV-2 genomic surveillance.

Genomic sequencing is essential to track the evolution and spread of SARS-CoV-2, optimize molecular tests, treatments, vaccines, and guide public health responses. To investigate the global SARS-CoV-2 genomic surveillance, we used sequences shared via GISAID to estimate the impact of sequencing intensity and turnaround times on variant detection in 189 countries. In the first two years of the pandemic, 78% of high-income countries sequenced >0.5% of their COVID-19 cases, while 42% of low- and middle-income countries reached that mark. Around 25% of the genomes from high income countries were submitted within 21 days, a pattern observed in 5% of the genomes from low- and middle-income countries. We found that sequencing around 0.5% of the cases, with a turnaround time <21 days, could provide a benchmark for SARS-CoV-2 genomic surveillance. Socioeconomic inequalities undermine the global pandemic preparedness, and efforts must be made to support low- and middle-income countries improve their local sequencing capacity.

Targeted genomic sequencing with probe capture for discovery and surveillance of coronaviruses in bats.

Public health emergencies like SARS, MERS, and COVID-19 have prioritized surveillance of zoonotic coronaviruses, resulting in extensive genomic characterization of coronavirus diversity in bats. Sequencing viral genomes directly from animal specimens remains a laboratory challenge, however, and most bat coronaviruses have been characterized solely by PCR amplification of small regions from the best-conserved gene. This has resulted in limited phylogenetic resolution and left viral genetic factors relevant to threat assessment undescribed. In this study, we evaluated whether a technique called hybridization probe capture can achieve more extensive genome recovery from surveillance specimens. Using a custom panel of 20,000 probes, we captured and sequenced coronavirus genomic material in 21 swab specimens collected from bats in the Democratic Republic of the Congo. For 15 of these specimens, probe capture recovered more genome sequence than had been previously generated with standard amplicon sequencing protocols, providing a median 6.1-fold improvement (ranging up to 69.1-fold). Probe capture data also identified five novel alpha- and betacoronaviruses in these specimens, and their full genomes were recovered with additional deep sequencing. Based on these experiences, we discuss how probe capture could be effectively operationalized alongside other sequencing technologies for high-throughput, genomics-based discovery and surveillance of bat coronaviruses.

C-to-U RNA deamination is the driving force accelerating SARS-CoV-2 evolution.

Understanding the molecular mechanism underlying the rampant mutation of SARS-CoV-2 would help us control the COVID-19 pandemic. The APOBEC-mediated C-to-U deamination is a major mutation type in the SARS-CoV-2 genome. However, it is unclear whether the novel mutation rate u is higher for C-to-U than for other mutation types, and what the detailed driving force is. By analyzing the time course SARS-CoV-2 global population data, we found that C-to-U has the highest novel mutation rate u among all mutation types and that this u is still increasing with time (du/dt > 0). Novel C-to-U events, rather than other mutation types, have a preference over particular genomic regions. A less local RNA structure is correlated with a high novel C-to-U mutation rate. A cascade model nicely explains the du/dt > 0 for C-to-U deamination. In SARS-CoV-2, the RNA structure serves as the molecular basis of the extremely high and continuously accelerating C-to-U deamination rate. This mechanism is the driving force of the mutation, adaptation, and evolution of SARS-CoV-2. Our findings help us understand the dynamic evolution of the virus mutation rate.

Insertion/deletion hotspots in the Nsp2, Nsp3, S1, and ORF8 genes of SARS-related coronaviruses.

The genome of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) contains many insertions/deletions (indels) from the genomes of other SARS-related coronaviruses. Some of the identified indels have recently reported to involve relatively long segments of 10-300 consecutive bases and with diverse RNA sequences around gaps between virus species, both of which are different characteristics from the classical shorter in-frame indels. These non-classical complex indels have been identified in non-structural protein 3 (Nsp3), the S1 domain of the spike (S), and open reading frame 8 (ORF8). To determine whether the occurrence of these non-classical indels in specific genomic regions is ubiquitous among broad species of SARS-related coronaviruses in different animal hosts, the present study compared SARS-related coronaviruses from humans (SARS-CoV and SARS-CoV-2), bats (RaTG13 and Rc-o319), and pangolins (GX-P4L), by performing multiple sequence alignment. As a result, indel hotspots with diverse RNA sequences of different lengths between the viruses were confirmed in the Nsp2 gene (approximately 2500-2600 base positions in the overall 29,900 bases), Nsp3 gene (approximately 3000-3300 and 3800-3900 base positions), N-terminal domain of the spike protein (21,500-22,500 base positions), and ORF8 gene (27,800-28,200 base positions). Abnormally high rate of point mutations and complex indels in these regions suggest that the occurrence of mutations in these hotspots may be selectively neutral or even benefit the survival of the viruses. The presence of such indel hotspots has not been reported in different human SARS-CoV-2 strains in the last 2 years, suggesting a lower rate of indels in human SARS-CoV-2. Future studies to elucidate the mechanisms enabling the frequent development of long and complex indels in specific genomic regions of SARS-related coronaviruses would offer deeper insights into the process of viral evolution.

A RT-qPCR system using a degenerate probe for specific identification and differentiation of SARS-CoV-2 Omicron (B.1.1.529) variants of concern.

Fast surveillance strategies are needed to control the spread of new emerging SARS-CoV-2 variants and gain time for evaluation of their pathogenic potential. This was essential for the Omicron variant (B.1.1.529) that replaced the Delta variant (B.1.617.2) and is currently the dominant SARS-CoV-2 variant circulating worldwide. RT-qPCR strategies complement whole genome sequencing, especially in resource lean countries, but mutations in the targeting primer and probe sequences of new emerging variants can lead to a failure of the existing RT-qPCRs. Here, we introduced an RT-qPCR platform for detecting the Delta- and the Omicron variant simultaneously using a degenerate probe targeting the key ΔH69/V70 mutation in the spike protein. By inclusion of the L452R mutation into the RT-qPCR platform, we could detect not only the Delta and the Omicron variants, but also the Omicron sub-lineages BA.1, BA.2 and BA.4/BA.5. The RT-qPCR platform was validated in small- and large-scale. It can easily be incorporated for continued monitoring of Omicron sub-lineages, and offers a fast adaption strategy of existing RT-qPCRs to detect new emerging SARS-CoV-2 variants using degenerate probes.